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Abstract 

We perform kinetic Monte Carlo simulations of flow-induced nucleation in polymer melts with 
an algorithm that is tractable even at low undercooling. The configuration of the non-crystallized 
chains under flow is computed with a recent non-linear tube model. Our simulations predict both 
enhanced nucleation and the growth of shish-like elongated nuclei for sufficiently fast flows. The 
simulations predict several experimental phenomena and theoretically justify a previously empirical 
result for the flow-enhanced nucleation rate. The simulations are highly pertinent to both the 
fundamental understanding and process modeling of flow-induced crystallization in polymer melts. 

PACS numbers: 64.60.qe, 64.70.km, 83.80.Sg 
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I. INTRODUCTION 



The nucleation of microscopic crystallites in polymer liquids is profoundly influenced by 
flow . This flow-induced crystallization (FIC), is a fascinating example of an externally 
driven, non-equilibrium phase transition, controlled by kinetics. Furthermore, FIC is ubiqui- 
tous in industrial processing of semi-crystalline polymers, the largest group of commercially 
useful polymers. A fundamental understanding of FIC promises extensive control of poly- 
mer solid state properties, as virtually every property of practical interest is determined by 
the crystal morphology. Flow can drastically enhance nucleation and trigger the formation 
of highly ahgned, e.ongated crystals, know, as shish kebabs Q. Recent experiments on 
entangled polymers, have studied, in detail, shish kebab formation [3J, |4], |5|] and the role of 
blend concentration {(J, molecular architecture [7] and molecular relaxation time f|. Often 
the most pronounced flow- induced effects occur near the melting point, where quiescent 
crystallization is immeasurably slow 

aa. 

The widely postulated mechanism for FIC states that flow forces the polymer chains into 
elongated configurations, which lowers the entropic penalty for crystallization [1]. How- 
ever, this hypothesis has yet to be developed into a quantitative molecular model. FIC is 
extremely sensitive to the flow-induced configurations of the non-crystalline chains, so an 
accurate molecular flow model is an essential prerequisite. Unfortunately, most polymer flow 
models predict only the macroscopic stress tensor and not the full molecular configuration. 
Alternatively, detailed simulations of polymer crystallization have provided much useful in- 



formation on the growth process [9|, [lfj, [ll|], yet simulating primary nucleation has proven 



difficult, especially at low undercooling, because of the extremely long nucleation times. At 
a much higher level of coarse-graining, models based on differential equations either assume 
an empirical dependence of the nucleation rate on the flow conditions 12], the stress tensor 

n n 

[131 ]. or the chain stretch [14|]; or assert that free energy changes under flow can be directly 
subtracted from the nucleation barrier 
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, Ji|. In either case the postulated FIC mechanism 
remains untested. An intermediate level of coarse-graining is required to surmount these 
difficulties. 



This letter presents coarse-grained kinetic Monte Carlo (MC) (17J simulations of 
anisotropic nucleation in flowing polymers. We compute chain configurations using a re- 
cent molecular flow model [l8j that reliably predicts both neutron scattering 19|, [20f| and 
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bulk stresses. Our simulations predict both enhanced nucleation and elongated shish nuclei. 
Kinetic MC has previously been used to model quiescent crystal growth in dilute polymers 



211 ] . However, it is particularly suited to nucleation and our algorithm is tractable even at 
low undercooling, providing an efficient and highly flexible framework to simulate general 
anisotropic nucleation under external fields. 

II. MODEL 

We compute the transient chain configuration of the uncrystallized chains under flow using 
;he GLaMM model jl8|. with finite chain extensibility included using Cohen's approximation 
22| | . The chains are divided into Z sub-chains, each corresponding to an entanglement 
segment of N e Kuhn steps of length b. We take N e = 100 throughout this Letter. One 
deterministic run of the model provides the end-to-end vector fj(t) = (r^), where the 
ensemble average is for sub-chains of type i. The data for an entire transient flow are used 
later in the nucleation simulations. All flow timescales are in units of the subchain Rouse 
time r e and r is normalized by 

Deformation of the amorphous chains has two effects on the nucleation kinetics: stretch- 
ing reduces the entropic penalty for crystallization; and monomer alignment modifies the 
probability of compatible alignment with the nucleus. The change in elastic free energy AF el 
for chains with ensemble average constraints f = (it), but locally at equilibrium, can be 



calculated by statistical mechanics [23j. Although an analytic calculation is not possible for 
finitely extensible chains, steep free energy gradients in highly stretched chains suppress fluc- 
tuations. Thus our numerical calculations for uniaxial deformations show that AF e '((rr)) 
can be accurately approximated by an expression that interpolates between Gaussian elas- 



ticity 23] for small Tr f and Cohen's 22jj approximation with r 2 = Tr f at high stretching, 



AF el = -Trf- -Tr lnf-iV e ln (l - — ) . (1) 



2 2 \ N t 

Similarly, numerical calculation of the monomer orientation distribution w(8) for chains with 
a constraint f are well approximated by using r 2 = Trf in the expression for w(8) derived 



from a direct constraint on r 



24]- 



w (0) = c-HVrn/N e ] h 



/Trf 



cos 8 
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where £ 1 is the inverse Langevin function and 9 is the angle between the monomer and 
the principle axis of f . Each subchain has an individual fj so is treated as a separate species 



with concentration 



1/Z. 



Shear flow 




FIG. 1: A sheared nucleus with protruding amorphous chains. 



Our coarse-grained simulations use the minimal nucleus description required for 
anisotropic nucleation. The nucleus comprises crystallized "monomers" or Kuhn steps, 
arranged in stems, with each stem formed from a single chain (see fig [I]). The total number 
of stems N s and the number of monomers in each stem is simulated. The arrangement of 
the monomers within the crystal is not resolved. The nucleus is assumed to be spheroidal 
with the polar radius L parallel to the stems. Assigning a crystalline volume of fejj to each 
monomer and normalizing all lengths by &o gives the equatorial radius W = \J N s /ir, the 
volume V = N T and thus the polar radius L = We also simulate the unit vector 

v parallel to the polar radius. As in classical nucleation theory the nucleus free energy 
comprises the free energy of crystallization, proportional to the nucleus volume, and a free 
energy penalty proportional to the surface area S. Thus the free energy in units of k^T is 
J-(Nt, N s ) = —€bN t + fisS, where e# and /x s are the coefficients of the volume and surface 
area free energies, respectively. 

We simulate two types of MC moves: addition of a new stem containing one monomer 
and lengthening of an existing stem by adding a new monomer, each with a corresponding 



reverse move. As in 



24] we assume that to attach, a monomer must be oriented within a 



solid tolerance angle Q of the nucleus orientation. For small Q, the fraction of monomers 
within this angle is w(8)Q, where 9 is the angle between the nucleus polar radius and the 
sub-chain principle strain axis. The stem attachment rate for species % is proportional to its 
melt concentration fa. In contrast, for stem lengthening only the next monomer along the 
chain forming the stem can crystallize. The concentration of this monomer at the nucleus 
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FIG. 2: (a) Transient and quasi-static nucleation rates for a Z = 25 monodisperse melt with 
6b = 1-9 and = 1.9 under start-up shear at -yT e = 0.1; inset contains GLaMM model predictions 
for shear stress (G e is the shear modulus), (b) Master-curve of nucleation rate against stretch ratio 
(e is the extension rate); inset shows the master-curves for varying en. (c) Steady-state nucleation 



251 ] compared with projected 



rate measurements against shear rate for an industrial iPP melt 
simulation results. 

surface where the lengthening event occurs is taken to be unity. Thus the attachment (+) 
and detachment (— ) move rates for stem addition (k st ) and stem lengthening (ki en ) are: 

k^ t = — <5W;min('l, e - ^ 8 "*) kZ = — min(T, e - ^^) 



St 



where r is the time for a monomer attachment attempt, Ui = Anwi{9i) is the fraction of 
correctly aligned monomers (normalized to unity in the quiescent limit), and the constant 
ln(f2/47r) has been added to ts- The free energy change of attaching a new stem of species 
% is 

= T{N T + 1, N s + 1) - F(N T , N s ) - ^Aif(f,,), 

where AFf is the flow-induced free energy change in sub-chain i. Similar calculations give 
the free energy changes for the other move types. A stem can only be detached if it contains 
a single monomer. 

The kinetic MC algorithm requires a sum over all possible move rates. The area available 
for stem addition moves is A(N S ), which is taken to be proportional to \^Ns~, and to give 
spherical nuclei in the quiescent limit. To obey detailed balance the rate of stem removal 
must be multiplied by A(N S — l)/N s , the probability of a given stem being at the nucleus 
surface. Each stem can lengthen or shorten from either the top of bottom. Thus the total 
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sum over all possible move rates is: 

Z , / „ r t \ n. 



tfxotai = A(N S ) <T(k%)i + A(N ^ 1] 

i=l 8 j=l 

3=1 

At each kinetic MC timestep one move is performed at random, with the selection probability 
weighted by the move rate. Time is then incremented by a stochastically determined interval 
At = — ln(y KTotah where ( is chosen uniformly on [0, 1] 13, \2l\ . 

After each MC step the nucleus orientation v is incremented over At by Brownian dy- 
namics. Flow rotates the spheroid nucleus through the Jeffery algorithm 26j|. The angular 
diffusion time scales with both the nucleus volume Nt and aspect ratio p following the 



expressions m 



261 ] . giving r rot = ar N T G(p), where G(p) depends only on the aspect ra- 



tio and we have introduced a dimensionless constant a, connecting r rot with the monomer 
attachment time tq. As the Jeffery algorithm is for Newtonian fluids we have neglected 
non- Newtonian effects here. After each time step all fj values are updated from the GLaMM 
model results, and AFf and u>i are recalculated. The parameter S = t £ /t sets the ratio 
of flow and monomer attachment timescales. We use S = 10 throughout this Letter, unless 
indicated otherwise, although similar results are obtained for any value of S > 1. 



III. RESULTS 



Each simulation evolves a single nucleus from a single monomer. The algorithm is espe- 
cially effective at low undercooling, as for small nuclei the rate sum is small, leading to large 
time steps. The quiescent free energy landscape can be calculated analytically from e# and 
Ps to give a dimensionless nucleation barrier Af* and a critical nucleus of n* monomers. The 
simulated nucleation time is the first time the polar and equatorial radii simultaneously 
exceed the critical radius r* = ^/3n* /An. Choosing a larger threshold size for nucleation 
has little effect on (tn). The results are accurately approximated by (rjy) = r exp(A/*). 
We obtained good statistics for barriers up to 25/cbT in ~ 50hrs on one 2GHz processor, 
giving a nucleation time of ~ lO 11 ^. The nucleation times are Poisson distributed, so the 
nucleation rate can be defined as N = 1/ (r N ). 

Under flow, the nucleation kinetics depend on the evolving chain deformation. We define 
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the instantaneous nucleation rate N(t) ~ T^m n ^ +A ^~"^~ A *^ where n{t) is the cumulative 
fraction of runs nucleated at time t. For shear rates 7 that are slow compared to the critical 
nucleus rotation time (7^7 <C 1) alignment effects can be ignored, equivalent to taking 
a = 0. The GLaMM model predictions for start-up of constant shear at 7r e = 0.1 of a 
Z = 25 melt and the resulting instantaneous nucleation rate are in fig[2](a). The results are 
independent of the ratio of flow and crystallization timescales for <S > 1. In fact, fixing the 
non- crystalline chain configuration to that corresponding to flow time t for the entire simula- 
tion and plotting the resulting quasi-static nucleation rate against t, reproduces the transient 
results (figEK). Thus nucleation is fully controlled by the instantaneous configuration of the 
surrounding chains, if r e > t . 

The stretch ratio A = A Yli=i V * s the dominant factor determining the nucleation 
rate, despite variations in flow rate, molecular weight and flow geometry (Fig[2](b)). This 
universal result is somewhat surprising as the distribution of stretch along the chain varies 
considerably with flow conditions, which may be expected to influence nucleation, espe- 
cially if molecular weight and flow geometry is varied. This result will be useful in deriving 
simple differential models of FIC (13I [yj] as the nucleation rate can be described the ex- 
pression N = iVo exp(?7(A 2 — 1)), where 77 is a fitting constant. Also in fig[2](b) is the curve 
N = N Q + /3(A 4 — 1), which empirically fits measured nucleation rates from flowing melts 
141 ] . The agreement with our simulation data for A < 3.5 theoretically justifies this empirical 
expression. In the inset to fig[2]^b) the sensitivity of nucleation to chain stretching increases 
with decreased undercooling (decreasing e#), as seen experimentally 22]. Fig[2]^c) directly 
compares our simulations with steady-state nucleation rate measurements on a polydisperse 
isotactic polypropylene sample (iPP) during shear 25j. The parameter determination, pro- 
jection of the simulation to large energy barriers and the approximation of the molecular 
weight distribution as a bimodal blend are detailed in appendix El The close agreement 
shows our model can quantitatively account for FIC measurements. Fig ]2]^a) also shows the 
transient nucleation rate when a = 5. This higher a value gives slower angular diffusion, 
meaning that flow aligns sub-critical nuclei, accelerating nucleation. Here, the quasi-static 
nucleation remains Poissonian and matches the transient rates, although iV now depends on 
both A and 7T™ D * t . Further increases of a give almost identical results. 

In experiments shish nuclei are especially prevalent in melts of short chains blended 

n it n □ n 

with a small amount of very long chains [J, [5|, la, 13, 18(. We simulated a melt of Z = 15 
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FIG. 3: Simulated nucleus aspect ratio at point of nucleation against shear rate for a 2% high 
molecular weight blend, (a = 5.0, /xg = 2.5, sheared for 120r e ); inset shows experimental orientation 
data on a bimodal blend sheared for 20 sec [g]. 

chains blended with 2wt% of Z = 52 chains, using a generalization of the GLaMM model 

n 

to bimodal blends [28]. High shear produces very elongated nuclei (fig [3]) from a purely 
kinetic mechanism. The shish widen by adding new stems using any monomer from the 
melt, whereas shish lengthen by adding monomers along an existing stem. Therefore the 
concentration of monomers from stretched segments at the growth surface is greater for 
lengthening than for widening, provided the nucleus contains a disproportionate number 
of stretched segments. Fast flow conditions are required for this disparity to overcome the 
significant surface area cost of elongated nuclei. Fig [3] (inset) shows experimental data for 
an orientational order parameter, the P2 orientation function, against shear rate for a 2% 
high molecular weight blend [8|, highlighting similarities with our predictions. Also in Fig [3] 
reduced undercooling, from reducing e#, increases the anisotropy, as seen experimentally 
j^ . When crystallization is less favorable, the disparity in the kinetics of stretched and 
unstretched chain segments increases. 



IV. DISCUSSION 



Our efficient kinetic Monte Carlo algorithm for flow-induced nucleation in polymer melts 
is tractable even at low undercooling. The flow-induced nucleation rate is a universal func- 
tion of the chain stretch ratio, independent of flow rate, molecular weight and flow geometry, 
but with decreasing undercooling causing increased flow sensitivity. This universal curve is 
very similar to an empirical relationship that accurately describes flow induced nucleation 
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rate measurements [14j . justifying this empirical result. Our successful quantitative compar- 



ison with nucleation measurements 25j establishes that changes in chain free energy under 
flow can describe flow enhanced nucleation in semi- crystalline polymers. In our bimodal 
blend simulations, a few percent of high molecular weight chains optimizes shish formation 
and the degree of anisotropy increases with shear rate and decreased undercooling, all of 
which are seen experimentally 

aaaas 

, |29J]. The simulation can readily be generalizec. 



to fully polydisperse melts, an essential step to model industrial polymer processing. The 
efficiency and flexibility of our algorithm makes it suitable to simulate general anisotropic 
nucleation under external fields. 

V. ACKNOWLEDGEMENTS 

We thank the EPSRC for funding (GR/T11807/01) and L Balzano, R Steenbakkers, G 
Peters, O Mykhaylyk, T Ryan and T McLeish for useful discussions. 

APPENDIX A: COMPARISON WITH NUCLEATION MEASUREMENTS. 
1. Introduction 

Here we quantitatively compare our simulation results with directly observed steady-state 
nucleation rates from a polymer melt under shear j^J . These measurements, by Coccorullo 
et. al, were made on a commercial grade isotactic polypropylene (iPP), known as iPP T30G. 
These data show an exponential dependence of the nucleation rate with shear rate, which 
our model correctly predicts. In order to quantitatively model this system, the required 
parameters fall into two classes: molecular relaxation times, which can be accurately calcu- 
lated from literature values; and crystallization parameters, for which order of magnitude 
estimates are provided by the literature, but whose specific values we determine from the 
nucleation data. The high polydispersity of iPP T30G is a further complication. As the full 
molecular weight distribution (MWD) for this resin was not reported, we choose a realistic 
MWD for highly polydisperse materials and demonstrate that this distribution is consistent 
with both the reported molecular weight averages and linear rheological measurements. For 
non-linear flow modeling, we approximate this polydisperse melt as a bimodal blend, since 
no polydisperse equivalent of the GLaMM mode is currently available. We note that no 



fundamental modification of our simulation algorithm would be required to take advantage 
of the results of a polydisperse tube model. In this data comparison we are able to correctly 
describe these steady state nucleation rate measurements against shear rate, using crystal- 
lization parameters that are consistent with estimates provided by literature measurements. 



2. Molecular relaxation times 



The relaxation times of iPP T30G can be determined by the tube model and confirmed 
against linear oscillatory shear measurements. Taking the tube model parameters for iPP 
from the literature {30J and shifting r e to the experimental temperature of 140°C using 
the time-temperature superposition parameters for iPP T30G [311 ]. gives r e = 90ns and 
M e = 4.4kg/mol. Here r e is the Rouse time of an entanglement segment and M e is the 
molecular weight between entanglements. The mass of a Kuhn step of iPP is 187.8g/mol 



321 ]. meaning the number of Kuhn steps per entanglement segment is iV e « 25. From these 



parameters the chain Rouse time can be calculated using tr = (M/M e ) 2 r e [331 ] . 

The form of the molecular weight distribution (MWD) for iPP T30G is essential to our 
modeling as small amounts of high molecular weight (HMW) material can control FIC (See, 



for example, ref [6|). Unfortunately, the full MWD of iPP T30G was not reported, with only 



the weight and number average molecular weights being given 



251 ] . Furthermore, standard 



distribution functions, such as the log-normal and generalized ex pon ential distributions, 



typically under-estimate the HMW tail of highly polydisperse melts [34J. Fortunately, linear 
rheological measurements are highly sensitive to this HMW tail at low frequencies and these 
data provide a clear upper-bound on the amount of HMW material. Thus we select a realistic 
form for the full MWD and show that this is consistent with both the reported complex 
viscosity and molecular weight measurements (M w = 481kg/mol and M n = 75kg/mol). We 
describe the MWD using a bimodal log-normal distribution, since the augmented tail of this 
distribution accounts more accurately for the HMW tail of highly polydisperse melts than 
standard monomodal distributions. This takes the form 

w B iend(m) = (1 - (p h )w LN (m, M wl , M nl ) + (p h w LN (m, M wh , M nh ) (Al) 



where <ph is volume fraction of the HMW peak, M w i/h and M n uh are the weight and number 
average molecular weights for the low and high molecular weight chains, respectively, and 
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wln is the monomodal log-normal distribution. Taking <ph = 0.01%, M w h = 3 x 10 5 kg/mol, 
M n h = 2 x 10 4 kg/mol, M w \ = 451kg/mol and M n i = 75kg/mol gives the correct overall 
weight and number average molecular weights and correctly predicts the complex viscosity 
measurements. These values were chosen as they correspond to the longest HMW tail 
that correctly predicts the complex viscosity measurements. Fig Ufa) shows the molecular 
weight distribution and fig H£b) compares the predicted complex viscosity resulting from 
this distribution with measurements on iPP T30G. Here the linear oscillatory shear response 
was computed using a polydisperse tube model for linear response, which combines the tube 
escape formula of Likhtman and McLeish 35[ with the Rubinstein and Colby algorithm 



for constraint release 
rheology, Rep Tate 



36| and is implemented within our in-house software for molecular 
41] . Using a monomodal log-normal distribution, with (pi = 0, produces 



slightly inferior agreement with the complex viscosity measurements and cannot account for 
the nucleation measurements. 



(a) 




M [kg/mol] M 




o Experimental Data 
r — Tube Theory Predictions 



C0[ 1/sec] 



FIG. 4: (a) Molecular weight distribution used to model iPP T30G calculated from eqn (jAljl . 

ngh molecular 
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The shaded area is the region of the molecular weight distribution defined as the 
weight tail in our modelling, (b) Linear rheological measurements on iPP T30G 
the predictions of polydisperse tube model using the molecular weight distribution discussed above 
and with G e = 1.1 x 10 6 Pa). 



As discussed above, because there is no suitable polydisperse tube model for non-linear 
flow, we approximate iPP T30G as a bimodal blend when modelling its non-linear response. 
FiglU^a) shows that there is no clear molecular weight at which to define the division between 
the two blend species. We choose the high molecular weight fraction to be all material with a 
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molecular weight above M* = 4200kg/mol, giving a HMW volume fraction of 1%. Taking a 
lower value of M* produces a Rouse time that is too short to account for the flow nucleation 
measurements, but we note that larger values of M* may be equally valid choices. We note 
also that this fraction contains discernible contributions from both of the monomodal terms 
in eqn flAlj) . We compute the expectation value of the Rouse time of the high molecular 
weight tail via, 

/ v J^T e (m/M e ) 2 w B iend(m)dm 

\ T R) = Foo p—r: , (A2) 

J M * w B iend{m)dm 

which gives tr = 66sec. All material below M* is defined as matrix material, and a similar 
calculation gives the average Rouse time of this matrix as 2.1ms. Thus, the matrix relaxes 
too quickly to be stretched by the flow rates in the study of Coccorullo et al.. Hence, we 
model the non-linear flow dynamics of this polydisperse material by approximating it as a 
bimodal blend with the two fractions having the Rouse times discussed above. 



3. Crystallization parameters 



To model iPP T30G we require its quiescent critical nucleus size and barrier height, at the 
experimental temperature. An estimate for the number of Kuhn steps in a critical nucleus 
can be found from literature data. The critical nucleus diameter has been shown to be 
of the order of the lamella thickness by transmission electron microscopy 
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for iPP 



and confirmed for a different polymer by atomic force microscopy 39j. For iPP the lamella 
thickness is ~ lOnm at 140°C. Using the iPP crystal density of 0.9g/cm 3 40j and the Kuhn 



step mass for iPP 32J and assuming a spherical critical nucleus gives the number of Kuhn 
steps in a critical nucleus as ~ 1000. 



The barrier height can be estimated from the nucleation data of Coccorullo et. al 25| (see 
fig [6]). Projecting these data to zero shear rate gives a quiescent homogeneous nucleation 
rate of 3.4 x 10 -11 //mi 3 /sec. This can be converted to a rate per Kuhn step using the melt 
density of iPP (0.85g/cm 3 ) and the iPP Kuhn step mass. Thus the density of Kuhn steps is 
Pk = 2.7 x 10 9 //im 3 so the quiescent nucleation rate is N = 1.2 x 10 _20 /sec per Kuhn step. 
Taking the crystallization timescale of a Kuhn step To to be of the order of the Kuhn step 
relaxation time tk = T e /N 2 = 0.144ns, means that NqTq ~ 10 -30 . Clearly nucleation in this 
system is extremely rare and such a separation of timescales cannot feasibly be simulated 
even by our highly efficient kMC algorithm. Instead we systematically project our results 
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to larger energy barriers to allow a comparison with these data. 



AI7k H 1 




FIG. 5: Projection of the simulated nucleation data to large energy barriers: (a) Simulation data 
for the average steady shear nucleation time against shear rate for a HMW volume fraction of 1% 
fitted with eqn (|A3j) [shapes are simulation data and lines are fits] ; (b) Variation of the sensitivity 
of the nucleation time to shear rate with quiescent barrier height; (c) Variation of tq with quiescent 
barrier height. 



4. Projection to large nucleation barriers 

To perform this projection we simulated the average nucleation time over a range of qui- 
escent nucleation barrier heights, varying from 10.6 to 26.3kBT, with the critical nucleus size 
held fixed at 1081 Kuhn steps. Nucleus alignment was neglected for all of these simulations 
(a = 0, throughout). The steady-state configuration for both blend components for various 
shear rates was calculated using a recent generalisation of the GLaMM model to bimodal 
The dimensionless steady-state nucleation time fjv = tn/tq against shear rate 



blends 



7 was simulated using the GLaMM model predictions, in steady state, for a range of shear 
rates and molecular weights. The results are independent of molecular weight when the 
shear rate is expressed as a Rouse Weissenberg number jtr, but do depend on the height 
of the quiescent nucleation barrier, A/*, as shown in figure 0(a). The simulated nucleation 
times vary exponentially with shear rate and can be described by the expression 

TnH) = fAf exp(-^7r fi ), (A3) 

where tno and £ are fitting parameters characterizing the decay of the nucleation time with 
shear rate, both of which depend on the barrier height Af*. The lines in fig 0(a) are the 
result of fitting eqn flA3j) to the simulation data. Both of the fitting parameters have a simple 
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dependence on A/*, as shown in figs E^b) and (c). Data for to the two largest quiescent 
barriers was used only to confirm the projection and was not involved in the fitting. This 
procedure allows the simulation results to be projected to nucleation barriers that are too 
high to be simulated. 



5. Data comparison 

We can use our projected results to model the data of Coccorullo et. al [25(. The value 
of the projected zero shear nucleation rate from these data, combined with the projection 
formula in fig[5](c) and the requirement that r ~ Tk, gives a quiescent nucleation barrier of 
A/* ~ TOkeT. The final value of /* is chosen to reproduce the slope of the experimental data 
in figure [6j Thus, using the Rouse time computed above and taking a value of A/* = 72keT 
gives £ = 0.33 from the projection formula in figure 0(b). This quiescent barrier specifies 
ttvo = 1-74 x 10 29 by the projection formula in figure [5](c). Thus taking To = 0.46ns produces 
agreement with the experimental data in fig [6] extrapolated to zero shear rate. We can also 
confirm that this value for r is ~ tr. With these parameters and using N(j) = 1/^(7), 



which is valid because all of the simulated nucleation processes are Poissonian, eqn flA3j) can 
be written as 



Pk 

Eqn (1A4I) produces the projected simulation data in fig El 



N(j) = _^exp(£7T fl ) = 3.4 x lO" 11 exp(21.8sec7)/sec//W. (A4) 
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? IG. 6: Comparison of measured steady shear nucleation rate against shear rate for an iPP resin 
with our projected simulation results. 

In summary, our model correctly predicts the exponential dependence of nucleation rate 
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on shear rate as observed by Coccorullo et al. 25j. Furthermore, with a characterization 
of the full molecular weight distribution and extrapolation of our simulation results to high 
nucleation barriers, we can make a direct quantitative comparison with these data. We use 
crystallization parameters that are estimated from the literature but whose final values are 
determined in response to the nucleation data. This results in very good agreement with 
directly measured nucleation rates under flow using parameters that are consistent with 
literature estimates. 

Several consequences arise from this data comparison. It demonstrates that changes in 
chain free energy upon stretching under flow are sufficient to quantitatively account for 
the degree of enhanced nucleation measured in flowing semi-crystalline polymers. It also 
indicates the strong need for comprehensive material characterization to accompany FIC 
measurements. In particular, an accurate characterization of the sample's high molecular 
weight tail, along with linear rheological modeling, will clearly define the full spectrum of 
molecular relaxation times. The issues of polydispersity could be solved either through 
FIC experiments on model polymers or through the emergence of a reliable non-linear flow 
model for polydisperse polymers. In principle, our model can predict the effect of varying 
molecular weight distribution, temperature and flow history, including transient flows, and 
a more comprehensive test of the model's predictions could be produced by systematic 
comparison to such a data set. The close agreement between theory and experimental data 
and our comparison provides a foundation on which further quantitative comparisons can 
be built. 
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